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Abstract. The ground state and thermodynamic properties of spin-1 and spin-| 
chains are investigated via exactly solved su(3) and su(4) models with physically 
motivated chemical potential terms. The analysis involves the Thermodynamic 
Bethe Ansatz and the High Temperature Expansion (HTE) methods. For the spin-1 
chain with large single-ion anisotropy, a gapped phase occurs which is significantly 
different from the valence-bond-solid Haldane phase. The theoretical curves for 
the magnetization, susceptibility and specific heat are favourably compared with 
experimental data for a number of spin-1 chain compounds. For the spin-| chain 
a degenerate gapped phase exists starting at zero external magnetic field. A middle 
magnetization plateau can be triggered by the single-ion anisotropy term. Overall, our 
results lend further weight to the applicability of integrable models to the physics of 
low-dimensional quantum spin systems. They also highlight the utility of the exact 
HTE method. 
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1. Introduction 

This paper is concerned with the study of two integrable quantum spin chains. By 
integrable we mean that the chains are exactly solvable in the non-diffractive PP Yang- 
Baxter sense jU 0] , with infinitely many conserved quantities and an underlying Bethe 
Ansatz solution. The two particular models are variants of the spin-1 and spin-| chains 
based on the family of integrable su(n) permutator models [HE]- The two-dimensional 
classical lattice model counterparts are the An-i models [HI Ej, aka the Perk-Schultz 
model [S] in a special case. 

Over the past few decades considerable effort has been invested in the analysis 
of exactly solved models in statistical mechanics OE]. Of particular relevance here 
is the calculation of the temperature-dependent properties [TT]. This began with the 
Thermodynamic Bethe Ansatz (TBA) approach ^21 1^ and has most recently evolved 
into the exact High Temperature Expansion (HTE) method ^lE]. The key ingredients 
of the HTE method are the Quantum Transfer Matrix (QTM) IT7j and the T- 
system ^H], from which one derives nonlinear integral equations which can be solved 
in an exact perturbative fashion. To date this approach has been applied to the 
Heisenberg model ^^j, the su{n) models ^2], the higher spin Heisenberg model fH] and 
to integrable quantum spin ladders [201 121] ■ In particular, it has been demonstrated 
that integrable models can be used to study real ladder compounds 120] • On the one 
hand, the TBA method was seen to be most convenient for predicting critical fields 
|221 12ni I21| of relevance to the quantum phase transitions in ladder compounds. On 
the other hand, the HTE method gave the temperature-dependent free energy, from 
which physical properties such as the magnetic susceptibility and the specific heat were 
derived. Excellent agreement between the theoretical and experimental results for the 
strong coupling ladder compounds has been found pUj . 

Here we apply the TBA and HTE methods to calculate the thermal and magnetic 
properties of spin-1 and spin-| chains.:}: 

1.1. Physical motivation 

Low-dimensional spin systems continue to reveal rich and novel quantum magnetic 
effects, such as fractional magnetization plateaux [221 12Z1 and spin-Peierls transitions 
[2H] • The effects of inter- and intra-chain interactions, single-ion anisotropy, biquadratic 
exchange interactions and the magnetic field all provide competing contributions in the 
low-temperature physics. As an example consider the spin-S* Heisenberg chain if 
25* is odd, the next-nearest-neighbour exchange interaction and anisotropy may drive 
the chain into either a non-magnetic singlet ground state (spin-Peierls transition) |2H| or 
a magnetic gapped ground state jHOI- Conversely, if 25* is even the Haldane gap may be 
closed in the presence of additional biquadratic terms or in-plane anisotropy. In general, 
the Haldane phase exists only in 25" even Heisenberg chains with antiferromagnetic 

X A brief account of our results can be found in Ref . [^ . 
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coupling, no gap appears in ferromagnetic Heisenberg chains. Accordingly one expects 
a rich quantum phase diagram to emerge for low-dimensional spin chains. More 
specifically, the intra- and inter-chain exchange interactions may result in either a 
valence-bond-solid Haldane phase or a dimerized phase jSI], whereas large single- ion 
anisotropy may trigger a field-induced gapped phase [HUl IS2] significantly different from 
the standard Haldane phase. 

The spin-1 Heisenberg chain has been extensively studied in the context of Haldane 
gapped materials [SSI OH IHSl IHEl Ell- Recently an anomaly in the susceptibility has 
been observed in the spin-1 compounds LiVGe206 |SH1- Such new materials have 
drawn attention from both theorists and experimentalists jSHl HHl 1^ ■ The Heisenberg 
spin-1 chain with additional biquadratic interaction has also been studied with regard 
to the Haldane phase [HI |l2]. In particular, the susceptibility of LiVGe206 can be 
quantitatively explained with the spin-1 bilinear-biquadratic chain [SHIIS]- On the other 
hand, a field-theoretic approach to the Heisenberg spin-1 chain with single-ion anisotropy 
was suggested by Tsvehk [22] ■ His theoretical predictions are in good agreement with 
experimental results for the compound Ni(C2H8N2)N02C104. 

1.2. Application to real compounds 

The Haldane gapped phase is observed in some one- dimensional materials exhibiting 
a weak single-ion anisotropy D, such that D < J, where J is the nearest-neighbour 
exchange interaction. However, a large anisotropy {D > J) can also produce a gapped 
phase 1301, which is significantly different from the Haldane phase found in weakly 
anisotropic compounds. In this latter case, the Haldane nondegenerate ground state is 
a valence-bond solid state in which a single valence bond connects each neighbouring 
pair to form a singlet. An expected excitation from the valence-bond solid state arises 
from a configuration in which a nonmagnetic state S'j = at site i is substituted for a 
state with Si = 1. In this way a total spin-1 excitation creates an energy referred to 
as the energy gap. However, bond alternation may lead to a dimerized ground state in 
which double valence bonds form two singlets between pairs of neighboring spins. By 
breaking one dimer with a nonmagnetic state at site i, a total spin-1 excitation is also 
separated. 

On the other hand, the field-induced gapped phase in the spin-1 chain is caused 
by trivalent orbital splitting, i.e., the singlet and doublet orbitals are separated by a 
crystal-field splitting. The singlet can occupy all states such that the ground state hes 
in the nondegenerate gapped phase. The lowest excitation arises as the doublet becomes 
involved in the ground state. This gapped phase is observed in some Nickel salts with 
large zero-field splitting, such as NiSnClg ■ 6H2O g^, [Ni(C5H5NO)6] (0104)2 El and 
Ni(N03)2 ■ 6H2O [IHj. Recently, other spin-1 magnetic compounds: 

(i) Ni(C2H8N2)2Ni(CN)4 (abbreviated NENC)|inillZl HE], 

(ii) Ni(CiiHioN20)2Ni(CN)4 (abbreviated as NDPK) gHl, and 

(iii) Ni(CioH8N2)2Ni(CN)4-H20 (abbreviated NBYC) 
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have also been identified as effective Heisenberg magnetic chains. This class of 
compounds exhibits a nondegenerate ground state which can be separated from the 
lowest excitation. Theoretical studies of these compounds have relied on a molecular 
field approximation via the Van Vleck equation [SU]. To first order Van Vleck 
approximation the exchange interaction is neglected and an effective crystalline field 
is incorporated to fit the experimental data. Not surprisingly these approximations 
cause some discrepancies in fitting the data. Although there has been some theoretical 
interpretation, the nature of the quantum phase transitions as well as the general 
microscopic Hamiltonian remains to be fully clarified for these compounds. 

1.3. Results and outline of this paper 

In this paper we consider integrable one-dimensional spin-1 and spin-| models with 
large planar single-ion anisotropy. We find that for the integrable spin-1 chain the 
non-magnetic singlet and the magnetic doublet states can be separated from the lowest 
magnon excitation by an energy gap. A large single-ion anisotropy together with in- 
plane anisotropy may trigger a non-magnetic mid-plateau, leading to a "ferrimagnetic- 
antiferromagnetic" phase transition. Significantly, a different type of gapped phase 
compared to the standard Haldane phase is found. We examine some real spin-1 
compounds via TBA and HTE. Excellent agreement between our theoretical predictions 
and the experimental results for spin-1 compounds such as Nickel salts and magnetic 
chains (NENC, NBYC and NDPK) is found. Our results show that the strong single-ion 
anisotropy, which is induced by an orbital splitting, can dominate the behaviour of these 
compounds. 

For the integrable spin-| chain the magnetic gapped ground state can be separated 
from the lowest magnon excitation, in contrast to the non-magnetic Haldane gapped 
phase in the spin-1 chain. A magnetization plateau at M = | originates at zero magnetic 
field. Moreover, an appropriately chosen value of the anisotropy can open a one third 
magnetization plateau reminiscent of the mixed spin-(l,|) ladder model [23] ■ 

The paper is organized as follows. In section |21 we introduce the integrable su(3) 
spin-1 chain and briefly review its exact solution. Then in section 12.31 we discuss the 
ground state properties via the TBA. In section 12.41 we derive the high temperature 
magnetic properties via the HTE method. We then discuss the application of the TBA 
and the HTE methods to real spin-1 materials with single- ion and in-plane anisotropics 
in section 121 In particular, we examine the thermal and magnetic properties of the Nickel 
salt NiSnCle ■ 6H2O (section ISU) and the magnetic chains NENC (section |S21), NBYC 
(section l3.3p and NDPK f section l3.4j) by comparing our TBA and HTE results with the 
experimental data for the specific heat, magnetic susceptibiliy and magnetization, where 
available. In section 0] we apply the same approach to the su(4) case in the context of 
the integrable spin-| chain. Concluding remarks are given in section [31 
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2. The integrable spin-1 chain 

2.1. The Heisenberg spin-1 chain 

Initially, the Haldane gap was identified in some quasi one-dimensional spin-1 materials 
|34| ESI IHE] described by the Heisenberg Hamiltonian 

L L L 

3=1 3=1 j=l 

in a parallel external magnetic field H. Here Si denotes the spin-1 operator at site i and 
L is the number of sites. The constants J and D denote the exchange spin-spin coupling 
strength and single-ion anisotropy, respectively. The Bohr magneton is denoted by /i^ 
and g is the Lande factor. For the compound NENP [22] the energy gap A between the 
ground state and the first excited state is of magnitude A ^ 12 K with the coupling 
constants J = 48 K and D = 0.2 J. While for the compound NINO A ^ 10 - 15 
K with J = 52 K and D = 0.3J. The above Hamiltonian has been recognized as a 
good model for Haldane-like compounds and has thus found widespread interest in the 
theoretical community [2niEIlE21- Unfortunately it does not appear to be integrable. 
However, if the spin-spin interaction term is changed, the model is integrable, and thus 
amenable to exact calculations using the sophisticated and well developed machinery 
of integrable models. We shall see below that under favourable conditions the modified 
Hamiltonian exhibits similar physical behaviour. 



2.2. The integrable su(3) spin-1 chain 
We consider the Hamiltonian 

H = J ^ iSj-Sj+i-\-(^Sj-Sj+^ j +H^i, 

cm. pot (2) 

i=i ^ ^ 
which only differs from the Heisenberg Hamiltonian (P) in the additional biquadratic 
interaction term.§ We discuss the on-site chemical potential terms Tichem.pot in detail 
further below. Leaving aside 7-^chom.pot, the model Q is the integrable sm(3) spin-1 
chain, which is well understood [SBIEII- It belongs to a family of solvable models with 
Lie algebra su{n) symmetry, which are generalized multi-state vertex models [71 133] . 
The spin-1 chain was originally considered by Uimin [3] in 1970 and later in the context 
of multi-state permutation operators by Sutherland 3]. A key ingredient is the identity 

Pi,, = Si ■ S, + {Si ■ S^,y - I, ® I,- , (3) 

§ There are in fact three integrable su(2)-invariant spin-1 Hamiltonians, each involving biquadratic 
interactions . Only the model considered here, with two-body interactions based on the permutation 
operator, appears to allow the desired chemical potential terms. 
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where I is the identity operator. We are thus considering the model 

L 

Ti. = J Pi,i + 1 + '^chcm.pot (4) 

i=l 

acting on V'^^ where is a three-dimensional vector space and Pi^i+i is the permutation 
operator, acting trivially on all sites, except sites i and i + 1, where 

Pi,i+i\ji, j2, ■ ■ ■ ■■■Jl) = J2, ■ ■ ■ , imJi,, • • • , Jl)- (5) 

We take the spin-spin exchange interaction J to be antiferromagnetic (J > 0). 

A crucial point in the following TBA and HTE analysis is the observation that 
the form of the permutation operator is not dependent on which basis spans V. 
Accordingly we always choose a local basis which diagonalizes the chemical potential 
terms. The chemical potentials characterize physical on-site interaction terms in the 
model Hamiltonian. Fortunately there is freedom to adjust them without losing 
integrability. 

To simplify the Bethe Ansatz, we apply periodic boundary conditions, i.e., we 
identify Pl,l+i = Pl,i and Vl+i = Vi- Following the standard approach this model can 
be solved by the nested Bethe Ansatz jH E] , which gives the energy eigenvalues 

^ = -'^Y.Vn)r. 7 + iVl/^l+iV2/^2 + iV3/i3 (6) 

,=1 {v) >f + ^ 

in terms of Bethe roots of flavour /colour Only the Bethe roots of first flavour 

fj^-* appear in the energy expression. The roots are determined by solutions to the Bethe 
equations 

Ml (1) (1) , . M2 (1) (2) 1- 
11 (1) (1) -11 (1) (2) , 1. ' 

1=1 - - 1 1=1 -^i +2^ 

Ml (2) (1) , 1- M2 (2) (2) , . 

TT % + 2' ^ TT n -^1+^ (J. 

11 (2) (1) 1- 11 (2) (2) .' ^'^ 

where j = 1, 2, . . . , Mi and A; = 1, 2, . . . , M2. 

In eq. © Ni denotes the occupation number of sites in state \i) for i = 1,2,3, 
where the basis |1),|2),|3) diagonalizes the chemical potential term under consideration 
with eigenvalues /ii, fi2, /^s- The Bethe quantum numbers Mi and M2, satisfying 
< M2 < Ml < L, are related to the number of sites in state \i) by 

A^i = L - Ml, 

A^2 = Ml - M2, (8) 
N3 = M2. 

The solution © in terms of the Bethe equations (|7j) is valid for any permutation of the 
basis order. We do not specify the basis here, rather we are free to choose the order of 
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the basis and we will exploit this in the TBA analysis of the ground state, as done in 
the TBA analysis of the integrable ladder models 



2.2.1. Chemical potential terms We choose the chemical potential terms which most 
closely model the physics of the compounds as used in the literature. The most general 
chemical potential terms exhibited in the real compounds are of the form 

L L L 

j=l j=l j=l 

L L L 

j=l j=l j=l 

The first term in each equation describes the effect of an external magnetic field H, 
which can be applied either parallel or perpendicular to the quantization axis given 
by the sample geometry. The in-plane anisotropy and planar anisotropy effects may 
be detected in magnetic properties of a single crystal. Actually, we can transform 
the transverse external magnetic field Hr^ in (fTUj) to the z direction via the following 
equivalent chemical potentials 

L L L 

j=l j=l j=l 

Therefore, the magnetic field perpendicular to the quantized axis provides a way to 
examine thermodjTiamic properties for the powdered samples via the chemical potentials 
flll|) . We will assume that D > and interpret it as a planar/single- ion anisotropy 
modelling the effect of crystal-field splitting. The third contribution describes an in- 
plane anisotropy effect which breaks the z"^ symmetry. We consider only the case E > 0. 

Although we do not need to consider the (complicated) eigenbasis of the chemical 
potential terms Q or (jlOp in performing the algebraic Bethe ansatz and applying the 
HTE method, the eigenbasis for different chemical potentials is helpful in analyzing 
the ground state properties of the Hamiltonian (0). The eigenbasis could for instance 
be presented as linear combinations of the spin projection onto the z-axis eigenstates 
15*^ = ±1,0) where the coefficients are dependent on H, D and E. The eigenbasis for 
the Hamiltonian (j2]) with different chemical potentials will be given in the following 
TBA analysis. The eigenvalues of the chemical potentials are given by 

/il = - ^/WTiil^^, (12a) 

fi2= -D, (126) 

/i3 = + if^sgHy, (12c) 

/il = - ^ - E + ^{D + Ef + {2fiBgHY) , {I2d) 

fi2= -E, (12e) 

1^3= - ^{d-E- ^{D + Ef + {2iiBgHY) (12/) 
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for and (jTU)), respectively. 

Note that we can adjust the chemical potential terms to model the physics of a 
specific compound without losing integrability, but we have no freedom in changing 
the spin-spin exchange interaction fixed by the permutation identity Q- In this way, 
the integrable sti(3) chain Q may describe some real compounds having either a spin- 
spin interaction close to the permutation interaction or exhibiting a strong single-ion 
anisotropy which effectively dominates the intrachain exchange interaction in the ground 
state. Fortunately, compounds are known where the relative strength J/D or J/E are 
very small. Thus the planar and in-plane anisotropics tend to dominate the physics. 

A large D anisotropy may trigger a gapped phase where the singlet occupies the 
whole ground state. In such a phase, the Bethe reference state becomes the true physical 
ground state, with Bethe numbers Mi and M2 in the Bethe equations ((Zj) equal to zero. 
The lowest excitation arises from a configuration in which the doublet is involved in the 
ground state. This excitation gives the energy gap, which can be caculated either from 
the energy (jH} in terms of the Bethe equations ((Zj) or from the TBA equations. We shall 
see that both the magnetic field and the in-plane anisotropy decrease the energy gap. 



2.3. TBA analysis: Ground state properties 

We consider the thermodynamic limit in which the Bethe equations ((Zj) have different 
types of solutions. The real roots, called magnons, form the ground state. Pairs 
of complex conjugated roots form bound states of two magnons. There are many 
kinds of bound states. The so-called string hypothesis was developed to classify such 
excitations, for example, for the one-dimensional boson gas with delta interaction |12j . 
the Heisenberg chain the Hubbard model (HZj, the supersymmetric t — J model [SHj 
and the Kondo model (23 • Applying the string hypothesis it is quite standard to derive 
the TBA equations 



-.00 00 
ln(l + r/W) = + 5^ A_ * ln(l + r^^"') - J] a_ * ln(l + r/^^)"'), 

m=l m=l 
^00 00 

ln(l + r^i^)) = -g'i^ + Y.A^^* ln(l + i^g^^) -J2<^nm* ln(l + (13) 

m=l m=l 

from ((Zj) in a lengthy procedure (TSj O ISHII221I23 1211 • We use the conventional notation, 
with 

1 Tl 

(^n{v) = 57— ^ (14) 

and the convolution denoted by *. The densities of roots and holes are pn\v) and 



for k = 1,2 with rjn^ = pn^^ (X) / pi!'\X) for n = 1, 2, 3, . . .. The driving terms 

(k) 

Qn depend on the choice of eigenbasis. We give them below for different configurations, 
i.e., we will permute the basis order to enable the simplest analysis of different regions 
in the parameter space (if, D, E). 
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Figure 1. Dressed energies for the su(3)-based TBA equations (|15ll . In each case 
the horizontal axis is the spectral parameter v with the origin {v — 0) at the minima. 
The vertical axis carries the dressed energies. Due to our basis reordering the upper 
graph (red) is always the current second excitation e*^^^ while the lower graph is the 
first excitation e^^-*. The Fermi energy for both flavours is always = 0, denoted by 
the dotted line. Shown are the four typical cases of solutions to (fT3|) : (A) no chemical 
potentials, (B) mixed state involving all three basis states, (C) gap in e*^^^ and (D) 
gapped ground state with energy gap A (see text in section 

We define the so-called dressed energies en^ via rjn^ := exp{en\X)/T) with k = 1,2. 
They play the role of excitation energies, measured from the Fermi energy level for each 
fiavour and satisfy Fermi statistics. In the limit T — 0, only the negative part of the 
dressed energies, e*^'-*" := min(e'-'\0), contributes to the ground state energy. Thus the 
TBA equations (fT^ reduce to two coupled nonlinear integral equations 

e(2)=^2 + ai*e('^--a2*e^'^-. (15) 

The solution of these equations provides a clear physical picture of the ground state 
properties. 

For typical solutions to (fTH|l see Figure To simplify the analysis of we 
separate the calculations for regions in parameter space where the basis order changes 
in accordance with the energy levels of the chemical potentials. The assignment \a) with 
a = 1, 2, 3 naturally depends on the parameters H, D and E, which change the energy 
level of the chemical potentials (as, e.g., in Figure |3)). The driving terms Qi are then 
positive or negative in defined regions [221 IIBI 121] allowing for a simple determination 
of critical fields and gaps. The band fillings are controlled by these driving terms. In 
the absence of symmetry-breaking potentials all bands are completely filled, i.e., their 
Fermi bundaries are their value at infinity, with e(oo) = ep. 
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2.3.1. D > 0, E > In this part of the TBA analysis we discuss the most general 
Hamiltonian with both anisotropies, D > and E > 0, with an external parallel 
magnetic field H. In later sections we specialize to cases without in-plane anisotropy 
to model certain physical compounds. The perpendicular magnetic field term ()10|) with 
eigenvalues (jl2(ij) -( []'2^ can be treated in the same manner. But the in-plane and planar 
anisotropies play different roles under perpendicular and parallel magnetic fields. 

Starting with Hamiltonian ^ without chemical potentials, i.e. with /i = 0, where 
the model has su{3) symmetry with the fundamental basis |1) = |2) = |0) and 

|3) = |— 1). The energy eigenvalue is threefold degenerate. The single-ion anisotropy 
term D{S^Y breaks this symmetry into a doublet and a singlet (sm(2) © u{l)). The 
fundamental basis is still an eigenbasis as the z direction is the axis of quantization. 
Turning on the external magnetic field, this symmetry is further broken into © 
u{l) © u{l). However, the in-plane anisotropy E breaks the symmetry even for 
H = 0. The energies split into three levels with respect to the anisotropy and magnetic 
field. The new basis, |1) = a_\-l) + |2) = |0) and |3) = a+|-l) + with 

a± = [fiB9H±^/{fiBgHy + E^]/E diagonalizes the potential terms © with eigenvalues 
()12aj) to ()12c|l . The basis order depends on the external magnetic field if we fix D and E. 
Note that still the original singlet state |0) remains unchanged. The two doublet states 
|±1) are rotated in the new basis. || To display the underlying physics in a clear way it 
is still convenient to use the singlet and doublet terminology, where the "doublet" now 
denotes the rotated former doublet states |±1), even if they are no longer degenerate. 

2.3.2. Case D > E In this section we discuss the configuration in which the single-ion 
anisotropy D is much stronger than the in-plane anisotropy E. The applied magnetic 
field is parallel as in (jHl). In this case the singlet |2) is initially energetically lower than 
the doublet, as per the left side of FigureEl The components |1) and |3) further split due 
to the Zeeman effect. In accordance with the energy levels of the chemical potentials, 
we change the basis order as the magnetic field increases. 

In the parameter region H < \/D'^ — E"^ / {hbQ) we have basis order |2), |3) with 
TBA driving terms 

91= --^ + D-^{iiB9Hf + E\ 

92 = 2,/{fiB9Hy + E^. (16) 
Whereas for H > \JD'^ — E"^ /{^b9) we have basis order |1), |2), |3) with driving terms 

9i= - -4vT -D + ^{liB9Hf + E\ 

92 = D + y/{fiB9Hy + E'. (17) 

In this way we see that the parameters D, J, E and H either raise or lower the Fermi 
surfaces of colours v'^^^ and v^'^\ A gapped phase occurs if all sites are occupied by the 

II Similar rotations have been considered by Huang and Affleck in a field theoretic study of the Haldane 
gap compound NENP (52) . 
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Figure 2. Chemical potential energy versus magnetic field for the parallel field case, 
equations H12a|l - H12cll . indicating level crossing. The left figure corresponds to equation 
()16|l and the right figure corresponds to equation H17(l . The solid line (black) denotes 
the singlet |2) with energy /i2, the dashed line (blue) denotes the doublet state |1) with 
energy /ii and the dot-dashed fine (red) denotes the doublet state |3) with energy ^3. 
In the left figure one level crossing occurs, marked by the dotted vertical line, while on 
the right side the basis order for the TBA analysis remains the same for all magnetic 
field. See text and also compare with Figure 01 



non-magnetic singlet in the ground state such that it takes a finite energy to excite the 
doublet. In the TBA picture this corresponds to the case where both excitation hands 
and are completely above the Fermi energy and therefore empty (see configuration 
D of Figure . From the TBA equations in the form (|TH|l this can conveniently be 
calculated by finding for which {H, D, E) the gap conditions 

ei{v = Q,H,D,E)>Q and e2(t; = 0, i/, D, E) > (18) 

hold. The restriction f = could be used because it can be concluded from the form of 
fll5|) that the absolute minima of the dressed energies occur at the roots f = 0. From 
(|TH|) and (fTKjl it can be concluded that a gapped phase lies in the region 

H<H,, = ^{D - 4 J)2 - E^fiBg. (19) 

Here Hd is the critical field indicating a quantum phase transition from the 
nondegenerate singlet to the Luttinger liquid phase. The energy gap is given by 

A = D - 4 J - TO^I^^^TrKE^. (20) 

This gap is clearly weakened by the in-plane anisotropy as well as the magnetic field. The 
necessary condition for this gap to exist can be read off from ()2()j) to be D > AJ + E. In 
the dressed energy picture this manifests itself by the point at which the lowest dressed 
energy e'-^^ just tips the Fermi energy = at its minimum, i.e. the transition D — > C 
in Figure HJ Note that due to our basis- reordering procedure the dressed energies always 
satisfy e^'^\v) > e^^^{v). Thus once H > Hd the component |1) with eigenvalue ()12ajl 
becomes involved in the ground state. At the critical point Hd the phase transition is not 
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of Pokrovsky-Talapov type, due to the mixture of the former doublet states \S^ = ±1) 
in the state |1) with a field-dependent magnetic moment. 

The magnetization increases as the magnetic field increases while more of 
component |1) becomes involved in the ground state. As usually defined in the TBA 
approach, the magnetization of an eigenstate of the system is 

3 

M{H) = J2N\,^rn\,^, (21) 

where the numbers N\j) denote the total number of state \j) with magnetic moment 

^b) = -^/^^■(^)- (22) 
Here the eigenvalues of the chemical potentials. The occupation numbers 

N\j),j = 1,2,3 are related to the Bethe quantum numbers Mi and M2 via (jH)), where 
Ml and M2 are given by 

Mi= [ p^^\v)dv and M2 = [ p^^\v)dv. (23) 



From the driving terms in ()17|) we see that once the field increases beyond the 
second critical point 

^ + (24) 

the ground state consists entirely of the |1) state. In this special case Mi = M2 = holds, 
and therefore the reference state |1) becomes a true physical state with (normalized) 
magnetization 

Mh>hM) = , ^""^^ . (25) 

We thus see that the magnetization in the pure |1) state gradually approaches the 
saturation magnetization M^ = 1 as the contribution to the magnetization from the 
original doublet state \S^ = —1) in the mixed state |1) tends to zero, whereas the 
state \S^ = —1) in the mixed state |3) becomes dominated. The probabilities of the 
components |1) and |2) are equal at the inflection point H = \JD'^ — E'^ / hbQ, where the 
magnetization Mjp = ^a/I — (E/D)'^. Actually, the magnetization between the critical 
fields can be obtained by numerically solving the TBA equations (I15|l . This reduces 
to solving one integral equation instead of two coupled integral equations because the 
second level dressed energy is gapful in many cases (i.e. e^'^\v) > 0). This numerically 
convenient case occurs when the general part B of Figure Q reduces to part C. Indeed, in 
Figure these novel phase transitions are observed in the magnetization derived from 
the TBA equations (fTH|) and HTE free energy given by (0H). As the spin-spin exchange 
interaction J decreases, the magnetization increases steeply in the vicinity of the critical 
point Hci- If J = 0, i.e., for the case of independent spins, the critical points Hd and 
Hc2 merge into one point, at which a discontinuity in the magnetization of the ground 
state occurs. This is in agreement with the trivial exact solution of a single-spin system 
with energy levels ()12aj) - p2c|l . 
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Magnetic Field H(Tesla) 



Figure 3. Magnetization versus magnetic field H in units of saturation magnetization 
for the Hamiltonian © with J = 0.2 K, Z? = 6 K, = 3.5 K and g = 2.0 with parallel 
magnetic field. The solid and dashed lines show the magnetization derived from TBA 
and HTE at T = and T — 0.3 K, respectively. The TBA magnetization curve 
indicates quantum phase transitions in the vicinity of Hd and Hc2 differing from the 
square root field-dependent critical behaviour. Their middle point has magnetization 
A/ip w 0.406 rather than A/jp = 0.5 in the case E = Q. Overall the TBA calculation is 
consistent with the HTE result. 



2.3.3. Case D < AJ + E If the single-ion anisotropy D lies in the region E — AJ < D < 
E + 4 J, the energy levels of the chemical potentials © are as shown on the right side 
of Figure 13 The gapped phase does not exist in this regime. There are two competing 
components in the ground state, namely |1) and |2). In this case the basis order is 
1 2), 1 3) and the driving terms 

91= ~ -D + ^E^ + {iiB9Hf 

92 = D + ^E^ + {fiB9Hy (26) 

apply throughout the whole range of the magnetic field. For very large magnetic field, 
i.e., H > He, where He = \/ {AJ + DY — E"^ / ^b9i the ground state consists purely of 
the doublet state |1). Again, in the dressed energy picture, H^ is the point at which 
the first excitation band, here given by the singlet state |2), touches the Fermi energy 
ejT' = at its minimum, as per the transition D — C in Figure H To the right, i.e., for 
higher magnetic fields, both excitation bands e*^^-* and e'-^'' are completely gapped (part 
D of Figure Q). However, in the case D < E — AJ, all three states are involved in the 
ground state in a gapless phase. 
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2.3.4- Perpendicular magnetic field If a perpendicular magnetic field is applied to 
Hamiltonian (j2)), similar quantum phase transitions might occur as long as the in-plane 
anisotropy is much stronger than the planar anisotropy. Namely, if > D + AJ, a 
gapped phase lies in the region H < Hd, where the critical field 



The magnetization increases almost linearly as the magnetic field increases above Hd- 
A quantum phase transition occurs at the critical point Hc2, where 

i/. = ^M±M(5±iZES. (28) 

The magnetization tends to saturation as the field becomes stronger than Hc2. However, 
if E < D + 4 J, there is no gapped phase even in the absence of the magnetic field. 



2.3.5. Special case 1: Absence of in-plane anisotropy (E = 0) Now consider 
Hamiltonian (j21) with chemical potentials Q for the special case E = 0. This will be 
used to analyze the compound NENC in section which is believed to have negligible 
in-plane anisotropy. The fundamental basis, |1) = |2) = |0) and |3) = |— 1), 

diagonalizes the chemical potential terms with eigenvalues 

/i2 = -D, 

/i3 = fJ'SgH. (29) 

Taking ii^ = in the driving terms in the different regions, we find that the ground 
state in the zero temperature limit has a gap if the single-ion anisotropy D > 4J. The 
singlet 1 2) ground state is separated from the lowest spin excitation by an energy gap 
A = D — 4J. This energy gap is decreased by the external magnetic field H. The 
pure singlet ground state breaks down at the critical point Hd = {D — AJ)/ ^Bgs- The 
magnetization almost linearly increases with increasing magnetic field due to the magnon 
excitation. We discuss this further via the HTE approach in section 12.41 The ground 
state is fully polarized once the magnetic field increases beyond the second critical point 
Hc2 = (-D + 4:J)/{fiBg), i.e., in the M = Ms saturation plateau region. A gapped phase 
can only exist for anisotropy values satisfying the 'strong anisotropy' condition D > 4 J. 
As shown in |22j the magnetization in the vicinity of the critical fields H^i and ifc2 
depends on the square root of the field, indicating a Pokrovsky-Talapov type transition. 
In this regime, the anisotropy effects overwhelm the biquadratic spin-spin interaction 
contribution and open a gapped phase in the ground state. 



2.3.6. Special case 2: Trivalent orbital splitting Here we consider a further 
specialization of the previous case Q which can be mapped onto it via the identification 

{H',D',E') = {H + 2Ji,D ~2Ji,0), (30) 
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Figure 4. Energy levels versus magnetic field for eigenvalues %\2\ of the operator 
(pTT|l . The solid line (black) denotes the singlet |0), the dashed line (blue) denotes the 
doublet 1+1) and the dash-dotted line (red) denotes the other doublet state |— 1). The 
vertical dotted lines and the roman numerals mark the change of basis order in the 
TBA analysis (see text). Shown are the three typical choices of coupling constants. 
For their connection to energy gaps and magnetization plateau see the text and also 
Figure El and Figure El 



where the primed constants are those apearing in the general model 0. We have not 
found an experimental spin-compound realization of this particular Hamiltonian so far. 
Nevertheless, we give the TBA analysis here, as it has new and interesting physics in its 
own right, as it incorporates a physical mechanism to open more unusal gapped states 
and magnetization plateaux. 

If we consider the orbital degrees in the Hamiltonian, trivalent orbital splitting, as 
well as interatomic effects, may result in an alignment of the internal magnetic field to 
prefer one component of the spin doublet. To represent this physical effect we add a 
projection operator to the chemical potentials, so that 

L L L 

H.,cho..pot = D Y^iS^f -JiY. ^J^^^ - 1) - ^'^9H (31) 
j=i j=i j=i 

This additional term gives rise to a rich phase diagram by shifting the energy levels, as 
seen in the difference between Figure |2] and Figure 01 The on-site energy eigenvalues 
are now given by 

/^l+i) = -fJ'BgH, 

m = -D, (32) 
/i|_i) = fisgH - 2Ji, 

where the fundamental basis, |1) = |+1), |2) = |0) and |3) = |— 1), diagonalizes the 
chemical potential term. 

In general, the magnetization curves for this case can have up to 3 gapped phases, 
as shown in part (B) of Figure O Note that all magnetization plots in this figure have 
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been obtained by the HTE method for relatively low temperatures and are therefore 
only similar to the zero-temperature magnetization. Numerical solution of the TBA 
equations gives slightly 'sharper' contours. Variations in the parameter set (J, Ji, 
D) can open or close the different gapped phases. After calculating the whole phase 
diagram via the TBA, analogously as done in Refs [221 ESI Elj , we show the physically 
interesting choices of coupling constants, corresponding to the four pictures in Figure El 
Here we have restricted our analysis to small exchange coupling values J < D/4. We 
now consider the model in the region Ji < D < Ji + 4 J, corresponding to configuration 
A of Figure 0] Proceeding with the basis orders 

(|3), |2), |1)) ^ (|3), |2)) -> (|1), |3), |2)) ^ (|1), |2), |3)), (33) 

we analyze the ground state properties from the TBA equations fll5|l with different 
driving terms. The technique is similar to the previous section and we omit the details 
of the driving terms. It is interesting that a degenerate gapped phase starts from zero 
field. The ground state consists of the doublet state |3) with spin 5*^ = — 1. This state 
is separated from the lowest magnon excitation by the energy gap 

A = 2 Ji - D - 4 J. (34) 

The magnetic field lifts the energy level of the state |3) and the gap vanishes when the 
magnetic field exceeds the first critical field 

„>„^,^ (35) 

The singlet state and the other doublet state |1) subsequently become involved in the 
ground state. Thus the magnetization increases almost linearly as the magnetic field 
increases. It is worth noting that the singlet state is suppressed by the strong potential 
Ji. The ground state is fully-polarized beyond the second critical magnetic field 

H>H^2 = (36) 

This phase diagram is shown in part (A) of Figure 

In the region Ji + iJ < D < 2Ji — 4J, we find the situation as depicted in 
configuration (C) of Figure ^ where the model exhibits three different gapped phases 
due to existence of a strong single-ion anisotropy. Here the basis orders become 

(|3), |2), |1)) ^ (|2), |3), |1)) - (|2), |1), |3)) - (|1), |2), |3)). (37) 

The plateau with total spin = —1 starts at zero field. The energy level of this state 
is increased by the external magnetic field, so that the spin singlet completely displaces 
the doublet from the ground state when the field exceeds H > Hc2- This field- induced 
nonmagnetic mid-plateau opens with the anisotropy parameter D. However, it can be 
closed for certain large values of J. The critical fields predicted by the TBA analysis 
are given by 

_2Ji-D-U _ 2 Ji - D + 4J 

Hci — , -nc2 — , 

fJsg fJsg 
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(38) 



D-AJ _ D + 4J 

J^c3 — ) -nc4 — • 



Further physical interpretation of these critical fields is given in section EiU via the HTE 
approach (see part (B) in Figure E))- If the single- iron anisotropy becomes larger, such 
that 2 Ji —AJ<D< 2 Ji + 4 J, the gapped spin 5*^ = — 1 state is weakened by the singlet, 
whose energy level decreases. Thus the doublet magnetic ground state is gapless as long 
as the field is below Hc2- In this case the mid-plateau can still survive (see part (C) in 
Figure El). Except for the absence of the first critical field ifci? the other critical fields 
in ()38|1 remain valid for this case. However, if D is large enough a nondegenerate singlet 
state can start from zero field. This scenario is depicted in part (B) of Figure 0] In this 
case, the anisotropy Ji is suppressed. The critical fields are the same as in section 1^.3.51 

2.4- HTE analysis: properties at high temperatures 

We begin by briefiy reviewing the idea of the HTE method, which we use to predict 
physical properties of the spin chain compounds in the next section. 

The Bethe Ansatz as described in section 2.1 is based on a row-to- row transfer 
matrix, which has the disadvantage that its largest eigenvalue may be degenerate. 
Moreover, the TBA equations suffer from infinitely many nonlinear integral equations 
(NLIE). On the other hand, the Quantum Transfer Matrix (QTM) jlS] considers the 
transfer matrix on a square lattice with an inhomogeneity parameter value mat = —f3J/N 
{N is the Trotter number). The largest eigenvalue is non-degenerate |T7] and the free 
energy is determined only by the largest eigenvalue of the QTM. The advantage of the 
QTM formalism is that the problem of deriving the largest eigenvalue reduces to solving 
finitely many NLIE. The nested Bethe equations take a slightly different form than in 
the row-to-row case. 

For the model at hand, the eigenvalue of the QTM is given by 



Q2{v) 



Qi{v-li) Q2{v 



+ e-^^^<f>4v)Mv + i)%^r^- (39) 

^2(^^-1) 

Here we follow notation from Ref. ^Hj, with (f)±{v) = {v ± iun)^ and Qaiv) = 
rii^i '("^ ^ — The Bethe roots are solutions of the QTM Bethe 

equations 

N 

,(1) I lA ^ Ml ^(1) „,(1) , : M2 „,(1) „,(2) 1. 



Vj - IUn + _ M1-M2 T-T - + 1 T-r - 2^ 

7i) '■ rl ~^ ^ 11~(T) (T) (2) , 1.' 
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11 „,(2) , 1. 11 „,{2) „,(2) .• '^^^^ 



, T 

(2) / 

' + iMTv + 1 / r=T ^fc" - + 7=1 ^k' - ^r' - 1 

In the above j = 1, . . . , Mi and = 1, . . . , M2. 

The eigenvalue t'^\v) is embedded in a family of eigenvalues Tm\v) for related 
models via the fusion hierarchy fH]- These eigenvalues are connected to each other by 
the functional relation ()42|) known as the T-system. The indices of Tm'* denote that 
the auxiliary space has fusion type (a, m). It carries the m-fold symmetric tensor of 
the a-th fundamental representation of the su(3) Lie algebra. The relations can be 
conveniently denoted graphically by Young tableaux, where the simple case only deals 
with rectangular shaped blocks (see Refs [T3] and J^). The fused eigenvalues are 
explicitly given by 

a m 

^ = E n n ^ - -m-2j + 2k)), (41) 

dj,k i=l k=l 

where the summation is taken over the "Young tableaux" parameter dj^k ^ {1, 2, 3} such 
that dj^k < dj^i^k and dj^k < dj^k+ii which form an a x m rectangular tableaux. After a 
renormalization the functional relation takes the simple form 

TirH^ + \^)Ti{:Kv - ii) = fi::i,{v)fitM + f^^~'\v)f^^^'\v). (42) 

in the normalized functions Tm^ 

By exploiting knowledge of analyticity properties, and numerically backed 
assumptions about the largest eigenvalue and the locations of its zeroes, this T-system 
can be transformed into a set of coupled NLIE. These NLIE had first been discussed 
for the spin-i chain (s-u(2)-XXX model) in ^1]. The general su{n) case was derived by 
Tsuboi jlSj in full detail. For the integrable spin-1 chain (su(3)-case) considered here 
the coupled NLIE are 
J 1 

dyTr(i/+/?r-ii)Tr(y+/5f'-ii) 



T^iv) = exp 

r«(.) = g«+ 



2vri (^;_^_/5«)r<^)(y+/5f'-i) L -(3[ 



1 -i> —1 



Ti^\v) = exp ((/ii + fi2 + f^3)/T) = const. (43) 

Here the contours are counterclockwise closed loops around zero, omitting certain 
special points 

The first remarkable property of these new NLIE's is that there are only finitely 
many of them, as the ?7,-th eigenvalue in this hierarchy is constant in terms of the spectral 
parameter v, in contrast to the finite temperature TBA considered in section EHl The 
second useful property is the fact that these NLIE's can be solved by an expansion in 
the small parameter J/T, for which the coefficients can be obtained recursively. This 
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distinguishes the HTE method from other NLIE approaches, e.g., Ref. ^7j, where the 
NLIE has to be solved numerically. 

The result for the free energy of the sm(3) case (already stated in Ref. ^1^) is 

- ^/(T, H) = In Q« + Cl, (0 + Cl, ' + Cl, (44) 

where the first few coefficients are given by 

n(2) 

Li n — 



1,0 — ^ .,n2' 

Qf Qf Qf 

in terms of the su(3) Q-system 

We find that considering only terms up to third order is sufficient for the analysis of 
real compounds. This is different from other types of high temperature series expansion 
|60j . The HTE coefficients are functions of the chemical potentials /Xj, i.e., they depend 
on the external parameters such as the coupling strength and magnetic field. The 
thermodynamic behaviour can be directly obtained from the free energy (j44j) via the 
standard relations 

M = -(^/(r,if))^, x = -[-gr.f{T.H)\, C = -T[^f{T,H))^, (46) 

for the magnetization, magnetic susceptibility and magnetic specific heat. We apply 
these results to several spin compounds in section El Now the small parameter of the 
expansion is J/T, i.e., the expansion is expected to work well for spin chain materials 
with relatively weak intrachain coupling. 

Before proceding, we note that in the above theory, the largest QTM eigenvalue is 
only fixed up to a (positive) constant. From the derivation of the HTE it can be seen 
that the free energy expansion 



/WT,i/) = Co + Ci- + C2(^-J +C3^^) +... (47) 
is only fixed up to a multiplicative constant, i.e. 

/phys. units 7conv. 



7c 



C0 + C1- + C2I-I +C3I-) +... 



(4^ 



CONTENTS 



21 



Physically, this amounts to the fact that in the theory we don't specify the unit of 
energy, i.e. there is no analogue of a 'characteristic length'. Mathematically it can be 
traced back to the multiplicativity of the underlying i?-matrix of the integrable model. 
In practice this factor depends on the choice of units for the physical free energy. 

Another remaining problem is the apparent widespread use of non-standard non- 
Si units in the experimental community for measurements on real compounds which 
sometimes even depend on compound specifics such as molecular mass. In eq. ()48|) the 
unknown constant 7conv. is a prefactor, i.e., the same in all orders of the expansion. The 
expansion is generated for interacting spins from the Bethe Ansatz. For zero order, the 
non- interacting case J = 0, it reduces to the exact trivial solution for L independent 
spins (here, e.g., for a 2-state spin-1/2 system): 

/o-thorder = 7conv.Co = -/CfiTlu ( c'^s^ + c'^B^ \ . (49) 

The second equality only holds for the special point J = 0, but as 7conv. is a constant it 
can be extended in eq. ipHj) to all J > 0. 

Once the free energy is thus known in SI units, say [J], it is trivial to obtain all 
other properties dervied from it in SI units via 

d 

-^^phys. units /'phys. units- (^^) 

The important point here is that HTE can make predictions in physical units, even if 
the original T-system is only valid up to a constant. Thus all conversion factors can 
in theory be calculated using the additional information from the non-interacting spin 
system. In addition, for real compounds the parameters D, E and J are given in Kelvin. 
For convenience, we omit in this paper the Boltzmann constant /c^ in the conventional 
notations D/kB, E/kB, J/ks- 

We found different spin-1 compounds have the same conversion factor for the 
susceptibility. In this way if we normalise the magnetic properties, such as the 
magnetization (M/M^), the conversion constant does not play any role in fitting the 
experimental data. On the other hand, the conversion constant for the specific heat is 
material dependent due to the choice of units in a given experiment. For convenience 
we state these factors in the relevant figure caption each time they are used. 



2.4-1 ■ Special case: Trivalent orbital splitting In this section we consider the special 
Hamiltonian ()31|) . already studied via the TEA method. Diagonalizing the Hamiltonian 
we find the chemical potential terms ()H2j) 

yUl = -jJ^BgH, 

/i2 = -D, (51) 

yU3 = fJ'BgH - 2Ji. 

We plot the magnetization versus magnetic field and the susceptibility versus 
temperature for different configurations in Figure El and Figure IHl corresponding to 
the parameter sets found in the TEA analysis in section 12. HI Comparing the results 
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Figure 5. Magnetization versus magnetic field at T = 0.8K for the Hamiltonian (0) 
evaluated directly from the HTE. The structure constants are (A) D = IIK, Ji = lOK 
and J = 0.5K; (B) D = 19. 5K, Ji = 13.5K and J = 0.5K; (C) D = 18K, Ji = 9K and 
J = 0.5K; (D) D = IIK, Ji = OK and J = 0.5K. In each case = 0.672K/T and 
5 = 2.0. The critical fields estimated from the magnetization curve coincide with the 
TBA results. 

of both approaches (see Figure El for numerical TBA values) we find that the critical 
fields obtained by the TBA coincide well with the values from the HTE free energy 
dm), even at relatively low temperatures. For case (A) in Figure El and Figure IHl 
where Ji < D < Ji + 4J, a field-induced gap exists in the lowest magnon excitation. 
The ground state hes in a magnetic doublet state with an energy gap. Note the typical 
rounded peak in the susceptibility for the antiferromagnetic gapped spin chain. For case 
(B), further enhancement of the anisotropy D into the region Ji + AJ < D < 2Ji — 4 J 
opens a non-magnetic singlet mid-plateau. This leads to a spin-Peierls-like quantum 
phase transition. The magnetic ground state is gapped, as before. If the single-ion 
anisotropy D is large enough, the gapped ground state is weakened by the singlet 
potential. If 2Ji — AJ < D < 2Ji + 4J, the ground state shifts into a gapless phase, 
as shown in case (C). The round peak in the susceptibility curve, indicating an energy 
gap, now gets replaced by a steep increase for temperatures close to zero, indicating 
a gapless phase. In this case, the mid-plateau remains open if the anisotropy is large 
enough, i.e., if 2Ji — AJ < D < 2Ji + 4 J. The most common gapped phase is found in 
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Temperature T(K) Temperature T(K) 

Figure 6. Susceptibility versus temperature at _ff = T. The numerical values are the 
same as in the previous figure. The susceptibility curves in (A), (B) and (D) indicate 
the existence of energy gaps, while no significant gap can be observed in (C). 

case (D), where Ji = and the anisotropy is large, D > 4 J. A non-magnetic singlet 
occupies the whole ground state as long as the energy gap A = D — AJ remains closed 
by the external magnetic field. We are now ready to move onto the direct comparison 
with experimental results. 

3. Examination of real compounds 

The compounds to be examined via the integrable spin-1 model (j21) are listed in Table 1. 
We summarize the coupling constants obtained from the Heisenberg chain (HC) and our 
model (HTE and TEA). We shall see below that the thermal and magnetic properties 
derived from the HTE method are in excellent agreement with the experimental results. 

3.1. Nickel salts 

Certain Nickel salts can be described by the spin-1 Hamiltonian (j21) with chemical 
potentials and (fTIHl with E = 0. The chemical potential terms are explicitly given 
in for the parallel case i.e. 

yUl = -fJ'BgH, 
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Compound 


HC J 


HC D 


HC E 


HTE J 


HTE D 


HTE E 


TBAiJci 


TBAifc2 


NiSnCl 


0.02K 


0.64K 


OK 


0.02K 


0.64K 


OK 


0.38 T 


0.49 T 


NENC 


0.65K 


6.15K 


0.7K 


0.17K 


6.4K 


OK 


3.8 T 


4.7 T 


NBYC 


0.2K 


2.55K 


1.5K 


0.35K 


2.62K 


1.49K 


N/A 


2.7 T 


NDPK 


0.96K 


5.0K 


2.8K 


0.55K 


5.15K 


2.55K 


1.0 T 


4.66T 



Table 1. Estimates for the exchange coupUng J and anisotropy constants D and E 
for some real compounds. Here HC denotes the values obtained from numerical fits for 
the spin-1 Heisenberg chain [461 l47lE5ll49j and HTE denotes the values obtained using 
the integrable spin-1 model and the HTE method. Also shown are the corresponding 
TEA estimates for the critical magnetic fields. 



= -D, 

and for the perpendicular magnetic field case ^TU^ . 

1^2 = 0, 

/i3 = -^P - VD^ + i2fiBgHy). 

A large zero-field splitting effect was identified in some nickel salts, such as NiSnClg • 
6H2O (abbreviated as NiSnCI) [Ni(C5H5NO)6](C104)2 H and Ni(N03)2 ■ 6H2O 
Here the Ni^"*" ion causes a large splitting between the singlet and doublet in the 
spin triplet in an axially distorted crystalline field. For D > and no magnetic field, 
the triplet is split into a lower singlet and an excited doublet. An early theoretical 
prediction for this kind of compound was done as a first order approximation via Van 
Vleck's equation jHO]- In this treatment the exchange interaction is neglected in first 
order for the specific heat, magnetization and susceptibility. An effective field was 
introduced for the magnetization and susceptibility to fit the experimental data. In 
section 12.31 we have shown that large anisotropics D may drive the antiferromagnetic 
spin chain into a gapped phase, which is significantly different from the valence-bond- 
solid ground state. Examining the magnetic properties of the salt NiSnCle ■ 6H2O jl3j 
via the HTE method, we find that the measurements are best fitted by a model with 
zero-field splitting D ^ O.GAK and a much weaker spin exchange interaction J ^ 0.027^^. 
Actually even the zero order HTE with free energy ()44|) gives a good analytic result for 
this compound As expected, the zero order term of the HTE free energy is just the 
result for L non-interacting identical spins with three energy levels fj,i, 1x2 and /^s, i.e. 

/(if,T),_,,,, = -Tln(e-^^/^ + e"^^/^ + e"^^/^). (52) 

Comparing this with the zero order HTE term in (03)? i-^-? with \n.Q^l\ we see that 
fl52|) is identical to the first term in the Q-System ()44|) . 

In FigurelZlwe compare the magnetization obtained from the theoretical free energy 
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Figure 7. Comparison between theory and experiment j43j for the magnetization 
versus magnetic field of the compound NiSnCle -6112 0. The magnetic field is parallel to 
the trigonal axis. A parameter fit suggests the coupling constants D = 0.64 K, J = 0.02 
K, 5|| = 2.22 with hb = 0.672 K/T. The conversion constant is Mrte « 8.5Mexp 
(10^ cgs/mol). The agreement between the theoretical curves derived here from the 
HTE approach and the experimental data suggests that the compound might exhibit 
zero-field splitting with much stronger splitting strength D than the spin exchange 
interaction J. 



expression (jl^ with parameters 

D ^ {)MK, J ^ 0.02ir, g\\ = 2.22, 

at temperatures T = 1.349ir, T = 2.162/^ and T = A.ISK with the experimental 
data for NiSnCl6-6H20. They show excellent agreement. The experimental curves are 
for the case of the external magnetic field parallel to the trigonal axis given by the 
crystal sample. The magnetic properties are anisotropic with respect to the external 
field direction. From our TBA analysis, at zero temperature, the energy gap is A 0.55 
K for this salt. At very low temperatures the system mainly occupies the non-degenerate 
singlet ground state if the magnetic field is less than the critical field Hd- The length of 
the antiferromagnetic correlations is finite. For higher temperatures T > D, the singlet 
ground state no longer dominates, with excitations caused by spin flips. 

In some experiments the samples are not of macroscopic crystal size, thus it is not 
possible to align the material parallel to an external field direction. In such 'powder' 
samples we assume an isotropic distribution and use the heuristic formula (see pl21 in 
Ref. |50j) 

Xpowder ~ lX\\ + lX±, (53) 
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Figure 8. Comparison between theory and experiment 1481 for the susceptibility versus 
temperature of the compound NiSnCl6-6H20. A parameter fit suggests the coupling 
constants D = 0.64 K, J = 0.02 K, g\\ = 2.22 and g±_ = 2.20 with /xs = 0.672 K/T. 
The conversion constant is xhte ~ 0.8f23xEXP (cgs/mol). The agreement between 
the HTE curve and the experimental data confirm that the compound NiSnCl6-6H20 
exhibits zero-field splitting with a larger strength than the spin exchange interaction. 



where x\\ and x± refer to the model Hamiltonian Q with respective on-site interactions 
dni) and (fTn|) . This simple ansatz describes the experimental data quite well, as can be 
seen in Figure |S] where all three susceptibilities - parallel, perpendicular and powder - 
have been measured for this compound. The temperature dependence of the magnetic 
susceptibility for powdered samples and the HTE predictions using the above averaging 
formula are given in Figure |H1 The typical rounded peak in the low magnetic field 
susceptibility curve, characteristic of low-dimensional antiferromagnets, are predicted 
to occur around 0.41 K from the HTE relation x\\ = —^TfHTE{T,H) for the parallel 
magnetic field. Using the same physical model, i.e., with the same parameters as for 
the specific heat comparison, we find that the susceptibility x\\ is in good agreement 
with the experimental curve. If the magnetic field is perpendicular to the trigonal axis 
(respectively along the x or y axis ), the susceptibility is given by x± = ~^2fHTE{T, H) 
with the perpendicular magnetic field chemical potentials. We see that our theoretical 
curve, depicted by the dashed line in Figure |H1 gives a good fit to the experimental 
susceptibility curve with perpendicular field. It also suggests the same parameter 
constants, except g± = 2.20. 

In Figure El the specific heat curve for H = T shows that the HTE result agrees 
with the experimental data [IH] in the available temperature region. Note that we 
do not change the coupling constants to get a best fit for each physical property. 
Rather we assign fixed coupling constants to each compound and thus truly predict 
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Figure 9. Comparison between theory and experiment |43| for the magnetic specific 
heat versus temperature of the compound NiSnClg -61120. A parameter fit suggests 
the couphng constants D ^ 0.64 K, J = 0.02 K, 5 = 2.2 with = 0.672 K/T. 
The conversion constant is Chte ~ 1.920CEXp(cal/mol — K) 8Cexp (J/mol-K). 
The dashed hne represents the experimental measured values Cp including the lattice 
contribution. The circle curve is the corrected magnetic specific heat without a lattice 
contribution of aT^ from Cp, where a w f .056 x 10~^ cal/mol-K"*. 

the physical behaviour. The dashed hne denotes the raw experimental magnetic specific 
heat Cp while the circles denote the experimental data, already corrected for a lattice 
contribution aT^ with a 1.056 x 10~^ cal/mol In the absence of a magnetic 

field a rounded peak, indicating short range ordering, is predicted around 0.23 K from 
the HTE result. Below T ^ 0.23 K, there is an exponential decay due to an ordered 
phase. 

3.2. The compound NENC 

It is known that antiferromagnetic spin-1 chains jSHESllSn! with weak planar anisotropy 
can exhibit a non-magnetic Haldane phase. The large D gapped phase has been observed 
in the compounds NENC, NDPK and NBYC |46| I49j . In these compounds the in-plane 
anisotropy — breaks the symmetry and weakens the planar anisotropy effect. It 
was inferred from the experimental analysis in NENC 0Bj that the contribution E from 
the in-plane anisotropy is negligible in comparison with large D, where the Nickel(II) 
orbit along the c-axis forms a strong crystalline field. As a result this strong crystalline 
field dominates the low temperature physics. The antiferromagnetic intrachain exchange 
interaction further lowers the energy but its contribution to the ground state, as well as 
the low-lying excitations, is minimal. As a consequence, the model with appropriate 
single-ion anisotropy is expected to describe this compound rather well. The specific 



CONTENTS 



28 




0123456789 10 
Temperature (K) 



Figure 10. Comparison between theory and experiment ^Hl for the magnetic specific 
heat versus temperature of the compound Ni(C2H8)2Ni(CN)4 (NENC). The conversion 
constant is Chte ~ SCexp (J/mol-K). The sohd line denotes the specific heat 
evaluated directly from the free energy with the parameters J ~ 0.17 K, = 6.4 
K, £■ = with g = 2.24 and hb = 0.672 K/T. The low temperature specific heat 
curves are shown in the inset. For low temperatures the inclusion of in-plane rhombic 
anisotropy E = 0.7 K (dashed line) gives a better fit than without rhombic anisotropy 
(solid line). 

heat has been measured up to a temperature around 10 K jlH]- A typical round peak 
for short range ordering at T 2.4 K is observed (see Figure ITUj) . An exponential decay 
is detected for temperatures below approx. 2.4 K. The calculated HTE specific heat 
evaluated from the model (j21) with J = 0.17 K, D = 6.4 K (the solid line with E = 
in Figure El) is in excellent agreement with the experimental curve in the temperature 
region T > 0.8 K. Our analytic result for the specific heat gives a better fit with the 
experimental data than the perturbation theory prediction For low temperatures 
(below 0.8 K) paramagnetic impurities and a small rhombic distortion are the main 
reasons for the discrepancy. The inset of Figure ^1 shows that the addition of a small 
rhombic anisotropy E = 0.7 K (dashed line) gives a better fit at low temperatures than 
that with E = (solid line). However, it is negligible for high temperatures. 

As far as we know the susceptibility measurement for this compound was made on 
powdered samples in a small magnetic field of if = 0.1 mT. The electron spin resonance 
for this compound was studied in Ref . j3H] • In Figure ^2 we compare the experimental 
results with our theoretical predictions. The susceptibility of NENC was measured only 
in the temperature range 50 mK - 18 K. From the data provided in Ref. jlH] we could not 
accurately estimate constants for the Curie- Weiss term or the paramagnetic impurities. 
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Figure 11. Comparison between theory and experiment for the susceptibility 
versus temperature of the compound NENC. The theoretical curve (solid line) is 
obtained via the empirical relation XPowder = \x\\ + |X-L ^-nd by adding a Curie- Weiss 
contribution. The inset shows the comparison for magnetization versus magnetic field 
at T — 4.27 K. The theoretical curve follows from the empirical expression for the 
powder magnetization. A fit for the susceptibility and magnetization suggests the 
coupling constants J = 0.17 K, D = 6.4 K and i? = 0, with = 2.18 and g\\ = 2.24. 
The conversion constants are the same as in the previous figure Figure |H| 



The theoretical HTE curves for the susceptibihty with both parallel and perpendicular 
magnetic field are again evaluated from the free energy ()44|1 with chemical potential 
terms Q and (jTUl). We obtained the powder susceptibility from the heuristic average 

XPowder = \X\\ + (54) 

In this form it does not fit the experimental data very well at low temperatures, probably 
due to small contributions from a Curie- Weiss term and paramagnetic impurities. We 
have corrected the powder susceptibility XPowder with the Curie-Weiss contribution 
c/ (T — 9) to obtain a better agreement with the experimental curves. Here c ~ 0.045 
cm'^ K/mol and 9 ^ —0.9 K. For our spin chain model, the CW-corrected experimental 
data suggests the fitted coupling constants J = 0.17 K, D = 6.4 K and E = 0, with 
g± = 2.18 and g\\ = 2.24. In this case the TEA analysis predicts an energy gap 
A ^ 5.72 K for a parallel external magnetic field. The inset of Figure ^2 shows the 
magnetization of a powder sample at T = 4.27 K. The singlet is supressed by the 
temperature. The best agreement between the experimental data and our powdered 
magnetization, Mpowdor = + is found with the anisotropic Lande factors 

g± = 2.18 and g\\ = 2.24. Here My and M± denote the magnetization for the external 
field parallel and perpendicular to the axis of quantization. 
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3.3. The compound NBYC 

In this section we discuss the apphcation of our theoretical results to measurements done 
on the spin-1 compound HH] Ni(CioH8N2)2Ni(CN)4 ■ H2O (NBYC). For this compound 
it is generally believed that the in-plane anisotropy E and the large anisotropy D play a 
significant role and so we expect the integrable model Q with suitable on-site anisotropy 
interactions to be a good microscopic model for this material. 

The thermodynamic properties have recently been experimentally investigated in 
Ref . jlH] • Theoretical studies based on a strong-couphng expansion method |62| suggest 
that the anisotropy of this compound might lie close to the boundary between the 
observed Haldane and field- induced gapped phases [IHl. However, due to the validity of 
the strong-coupling expansion method there are apparent descrepancies in fitting 
the specific heat, susceptibility and magnetization if the rhombic anisotropy is large. 
Here we evaluate the relevant quantities for the powdered samples using the free energy 
expression pij) with parallel and perpendicular magnetic fields, in-plane anisotropy E 
and single-ion anisotropy D (see and (fTn|l and their eigenvalues ()12ajl -( [T2^ ) via the 
empirical formula XPowder = + 

In Figure ^1 we compare the results for the susceptibility of NBYC. A fit for the 
susceptibility suggests the coupling strengths D = 2.62 K, E = 1.49 K and J = 0.35 K, 
with g\i = g± = 2.05. The small discrepancy at low temperature is probably caused by a 
Curie- Weiss contribution. The inset of Figure IT^ shows plots of the magnetization of the 
powdered samples at T=5, 10, 20 K. Again our theoretical curves are evaluated from the 
empirical relation Mpo„der = |^|| + An overall agreement for the magnetization 

at different temperatures is consistent with the parameters used for the susceptibility. 
The singlet state is now supressed by both the in-plane rhombic anisotropy and the 
temperature. 

We turn now to the specific heat which has been measured up to a temperature of 
6 K |49j. The theoretical specific heat evaluated from the integrable model (j2I), with the 
same parameters given above, is in good agreement with the experimental data in the 
temperature region 0.5 K to 6 K, as can be seen in Figure El For lower temperatures 
the HTE does not give reliable results due to poor convergence. 

3.4. The compound NDPK 

Another compound considered in Ref. jlHj to be a spin-1 magnetic chain with planar and 
in-plane anisotropy is NDPK. As far as we are aware, there has been no comprehensive 
theoretical and experimental study of this compound. The specific heat of NDPK was 
measured only in the temperature range from 100 mK to 2.5 K. A numerical ESCA study 
of the specific heat indicated this compound lies in the large-D phase. Quantitative 
agreement with the experimental data suggested the values D = 5 K, J = 0.96 K and 
E = 2.8 K We examined the specific heat via the integrable model (j21) with the 
parameters D = 5.15 K, J = 0.55 K and E = 2.55 K. The specific heat (solid line ) 
evaluated from the HTE is in good agreement with experimental data in the temperature 
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Figure 12. Comparison between theory and experiment 49 for the susceptibihty 
versus temperatm'e of the spin-1 chain compound NBYC. The conversion constants 
are the same as for the above materials. The solid line is the susceptibility evaluated 
directly from the free energy of the spin-1 chain for the powdered samples with coupling 
constants J = 0.35 K, D = 2.62 K and = 1.49 K, with g\\ = g± = 2.05. The small 
discrepancy at low temperature might be attributed to a Cuire-Wess contribution. The 
inset shows the magnetization for the powder samples at T=5, 10, 20 K. 

region 0.9 K to 2.5 K (see Figure IT^ . For temperatures below 1 K, the HTE is not valid 
due to the small ratio of J/T. We believe that the integrable model will also describe 
other magnetic properties of NDPK. 

4. The integrable spin-| chain 
4.I. The exactly solved su{A) model 

It is generally believed that the ground state of the antiferromagnetic spin-S" Heisenberg 
chain remains in a gapless phase when 25* is an odd integer fl^ I31j . Nevertheless, an 
appropriate projection operator can open a "Haldane gapped phase" [lOj. Numerical 
study [HI] of the spin-| antiferromagnetic chain has verified the existence of a one third 
saturation magnetization plateau. In this section we consider an integrable spin-| chain 
with two different kinds of single-ion anisotropics, capable of triggering gapped and 
gapless phases. We investigate the thermodynamics of the model in terms of the HTE 
and TEA. The two Hamiltonians in the fundamental basis are given by 

L L L 

^ = E + ^ E - ^^9H (55) 
j=i j=i j=i 
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Figure 13. Comparison between tlieory and experiment ^HI for the magnetic 
specific heat versus temperature of the compound NBYC. The conversion constant 
is Chte ~ IOCexp (J/mol-K). The solid hne denotes the specific heat at _ff = 0.1 mT 
evaluated directly from the free energy with parameters D — 2.62 'K, E — 1.49 K and 
J = 0.35 K, with g = 2.05 and hb = 0.672 K/T. 



where 

0=0 Pt^Ol ' 

with Xa = \a{a + 1) — S{S + 1) We define model 1 with projection operator 
and model 2 with P~, where 

^ = (5T + i)(5I-|)(57 + i), 

Pi = iS]?. (57) 

In the above, Sj is the spin-| operator acting on site j, J is the exchange coupling 
and D is the single-ion anisotropy. The operator projects spin states onto the state 
with = i, which leads to a magnetic ground state. In general, one can incorporate 
different projection operators in the spin-S* chain associated with the integrable su{n) 
model with n = 2S + 1 56j such that multi-plateaux phases can occur. 
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Figure 14. Comparison between tlieory and experiment |46| for the magnetic 
specific heat versus temperature of the compound NDPK. The conversion constant 
is Chte ~ IOCexp (J/mol-K). The solid Hne denotes the specific heat at 7? = T 
evaluated directly from the free energy with the parameters D — 5.15 K, E — 2.55 K 
and J = 0.55 K. 



4-.2. HTE approach 

To investigate the thermodynamic properties we again apply the HTE scheme. The 
eigenvalue of the QTM (up to a constant) is [SH E| 



Qi{v - 



Qi{v - i^i)Q2{v - 1) 
Q2KV - i)Qz[v - fi) 



+ e-^-V_(^)04- + i)^7M^- (5 



2 

5; 



In the above equations, /ij, i = 1,2,3,4, are the eigenvalues of the on-site chemical 
potentials with regard to the fundamental basis \S^) 

|1) = ||), |2) = |i), |3) = |-1), |4) = |-|). (59) 

For model 1 with projection operator P+ the chemical potentials are 

/ii = - IfisgH, /i2 = -^nsgH - 2D, 

/^3 = IfJ'BgH, /i4 = jfisgH. (60) 
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For model 2 with projection operator P the chemical potentials are 
/ii = - IfiBgH, fi2 = -\iJiBgH - 2D, 

/W3 = IfJ'BgH - 2D, /i4 = IfisgH. (61) 
The HTE expansion of the free energy, up to third order, is given by 

- ffiT, H) = In g« + + Cl, ' + Cl, ^£ 

where now [1^1 



+ •■■.(62) 



.(1) _ o ^1 



(2) 
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Vi Vi 

Here the Q-system for su{A) is 

_ g-/3Mi-/3/i2 _|_ ^-f^fii-f^fis _|_ g-/3/ii-/3^4 _j_ g~/3M2~/3M3 _j_ g-/3At2-/3/i4 _j_ ^-I3ij.s-l3fi4 

_ Q-'l3fJ.i-f^fl2-f3fJ.3 _|_ g-/3Ml-/3M2-/3/i4 _|_ Q-f^fll-l3fl3-f^fl4 _|_ ^-f3fl2-f^fl3-l3fj.4 

g(4) ^ g-/3Mi-/3/^2-/3M3-/3A'4_ 

^.c?. Ti?^ analysis: ground state properties and magnetization plateau 



± 



First consider Hamiltonian with projection operator The anisotropy terms P 
and the external field can be diagonalized simultaneously in the fundamental basis (jSUj) . 
If D > 0, the projection operator energetically favours the component |2). Proceding 
in an analogous fashion to the spin-1 TBA analysis, we find that the magnetic ground 
state is separated from the lowest magnon excitation by the gap A = 2D — 4J in the 
zero temperature limit. A magnetization plateau with M = | starts from zero field. 
The gap is diminished by the magnetic field. If the magnetic field is larger than the 
critical field Hd = {2D — 4:J)/fiBg, a phase transition from a ferrimagnetic phase to 
an antiferromagnetic phase occurs. As the magnetic field increases, the magnetization 
almost linearly increases up to the saturation magnetization. The saturation critical 
field is Hc2 = {2D + AJ)/ fiBg. The magnetization and the susceptibility evaluated from 
the free energy ()62j) with the chemical potentials given by ()60|) are shown in parts (A) 
and (A') of Figure ITHl The magnetization curve indicates that the ground state remains 
gapped until the magnetic field exceeds the critical field Hd- A round peak is found in 
the susceptibility. This conclusion coincides with the suggestion by Masaki et al. (SUl 
that the gapped phase might exist in the 25" odd integer antiferromagnetic spin chains. 

Apart from this degenerate ground state, there also exists another gapped M = | 
magnetization plateau in the spin-| chain. If we consider the Hamiltonian (j55j) with 
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a large anisotropy term P~, the ground state lies in a gapless phase. The gapless 
ground state could be shifted into the gapped phase with a magnetization plateau 
M = The field- induced magnetic effects occur as follows: (i) the system undergoes a 
phase transition from an antiferromagnetic phase into a ferrimagnetic phase at the 
critical point Hd = AJ/hbQ] (ii) the M = \ magnetization plateau vanishes at 
Hc2 = (2D — 4:J)/ fiBg; (iii) the ground state is fully-polarized beyond the critical field 
Hc3 = {2D + AJ)/ fisg- The magnetization curve given in part (B) of Figure El indicates 
the existence of a mid-plateau, which is reminiscent of the mixed spin-(l, i) ladder |2S1- 
Remarkably, with regard to this mid-plateau, a cusp-like phase transition emerges in 
the susceptibility curve in (B') of Figure UHl This field-induced cusp is similar to the 
spin-Peierls transition in spin-i chain materials [2H1- We anticipate that this novel phase 
transition may exist in real spin-| chain compounds. 

5. Conclusion 

We have studied the thermal and magnetic properties of spin-1 and spin-| chains via the 
su{3) and su{A) integrable models. This has enabled some novel phase transitions and 
critical behaviour to be predicted via the TBA analysis and HTE method. In particular, 
the zero temperature TBA predictions and the temperature-dependent analytic result 
derived from the HTE provide the full ground state properties and thermodynamics 
of the underlying models. The integrable spin-1 model with appropriate chemical 
potentials gives a good description of some real spin-1 chain compounds with large 
single-ion anisotropy. The planar anisotropy on-site interaction leads to a gapped 
phase which is significantly different from the Haldane phase in which the spin-spin 
exchange interaction results in a valence bond solid ground state. We have compared 
the theoretical curves with measurements on spin-1 chain compounds such as the Nickel 
sah NiSnCle ■ 6H2O (section EH}, NENC (section EH, NBYC (section ESI, NDPK 
(section lH.4j) . In general, the agreement with the experimental data for the compounds 
under comparison is excellent. Moreover, a specific anisotropy in the spin-1 chain can 
trigger a mid-plateau which causes a spin-Peierls-like quantum phase transition. 

We have considered two types of on-site anisotropies for the integrable spin-| chain. 
In contrast to the non-magnetic gapped phase in the spin-1 chain, the magnetic gapped 
ground state can be separated from the lowest magnon excitation in the integrable spin- 
I chain. Counting the magneitzation plateau in units of the saturation value Mg = |, a 
one third magnetization plateau with M = | can start at zero magnetic field. Moreover, 
an appropriate single-ion anisotropy can open a one third mid-magnetization plateau 
as also exhibited in the mixed spin-(l, i) ladder model [22] ■ We anticipate that the one 
third magnetization plateau may be observed in some real spin-| compounds. 
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Figure 15. Theoretical curves for the spin-| chain model (|55|l with single-ion 
anisotropy. Parts (A) and (B) show the magnetization versus magnetic field at 
T = 0.8K. In part (A) D = 5 K and J = 0.5 K which suggests a magnetic gapped 
ground state for the Hamiltonian H55|l with projection operator P+. In part (B) D = 6 
K and J = 0.4 K which suggests a gapless ground state for the Hamiltonian H55(l 
with the single-ion anisotropy . We have set fis = 0.672 K/T and g = 2.20. The 
critical fields estimated from the magnetization curves coincide with the TBA results. 
Parts (A') and (B') show the susceptibility versus temperature at iJ T. The 
numerical values in (A') and (B') are the same as in parts (A) and (B), respectively. 
The susceptibility curves in (A') indicate the existence of an energy gap, whereas no 
significant gap can be observed in figure (B'), however, the mid-plateau suggests a 
spin-Peierls-like quantum phase transition. 
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